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Recently it has been shown that time optimal quantum computation is attained using the Cartan 
decomposition of a unitary matrix. We extend this approach by noting that the unitary group is 
compact. This allows us to reduce the execution time of a quantum algorithm Ua, further by adding 
an extra gate W to it. This gate W sends Uag to another algorithm WUaig which is executable in 
a shorter time than Uag. We call this technique warp-drive. Here we show both theoretically and 
experimentally that the warp-drive reduces the execution time of Grover’s algorithm implemented 
with a two-qubit NMR quantum computer. Warp-drive is potentially a powerful tool in accelerating 
algorithms and reducing errors in any realization of a quantum computer. 
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Quantum computing is an emerging discipline based on 
encoding information into a quantum-mechanical system 
ua. There a quantum algorithm is expressed in a form 
of a unitary matrix 


i rT 
Vas = Ton |= f HOO, (1) 


where 7 stands for time-ordered product, y(t) denotes 
collectively the control parameters of the Hamiltonian H 
at time t. This unitary matrix is often implemented in 
terms of so-called elementary gates, such as U(2) gates 
and CNOT gates a [al E It is possible, instead, to 
directly implement a given unitary matrix without de- 
composing it into these elementary gates. It is expected 
that this will reduce the execution time required in gen- 
eral. One method for direct implementation is to employ 
numerical optimization of the control parameters of the 
Hamiltonian 6, [71 [g]. The other method is to use the 
Cartan decomposition of the group SU(2”), to which an 
n-qubit matrix representation of a quantum algorithm 
belongs g id. The latter approach, which is adopted 
here, is successfully demonstrated recently using an NMR 
quantum computer, whose pseudopure state is generated 
by cyclic permutations of state population fia}. It should 
be noted that exact optimal implementation of a quan- 
tum algorithm has been achieved in holonomic quantum 
computation in an idealized case ug. 

A warp-drive is a fictitious gadget with which two re- 
mote points in space are connected Ep It is the purpose 
of the present Letter to demonstrate, both theoretically 
and experimentally, that a similar technique may be em- 
ployed to shorten the execution time of a quantum algo- 
rithm. A time optimal implementation of the quantum 


algorithm is equivalent to navigating along the time op- 
timal path from the identity operator to the point Uag 
in SU(2”) by tuning the parameters y(t) in the Hamil- 
tonian. We show below that an additional permutation 
matrix W of the basis vectors, when added after Uag, 
sends it to a point WUag near the identity matrix I 
so that it takes a shorter time to follow the time opti- 
mal path connecting I with WUag than I with Uag, see 
Fig. 1. Therefore we have “warp-driven” Uag to WUag 
by adding W. Although it might seem counterintuitive 
that the execution time is reduced by adding an extra 
gate, this is the case since the unitary group SU(2”) is 
compact. Let us consider navigating on a sphere $°, 
which is also a compact space. Suppose we leave from 


FIG. 1: Conceptual diagram showing the effect of a warp- 
drive gate W in a compact space SU(4). A matrix Uag is 
sent to WUalg which is reachable from the unit matrix J in a 
shorter execution time. The curves connecting these matrices 
represent time optimal paths. 


a point O on the equator toward west. Then we will 
eventually arrive at the antipodal point. The distance 
between O and us will be shorter and shorter as we fur- 
ther circumnavigate the sphere. This is what happens 
when an extra gate is added. 

In an NMR quantum computer, a one-qubit operation 
may be carried out in a short time on the order 10 ps 
while a two-qubit entangling operation takes time typ- 
ically ~ 10ms. Thus one-qubit operation time may be 
ignored in estimating the overall execution time. Let us 
consider a molecule with two heteronucleus spins for def- 
initeness, whose Hamiltonian, in the rotating frame with 
respective Larmor frequency, is 


H(y) = —w [cos ¢1(o2 ® 12/2) + sin di (ay ® I2/2)] 
—W12 [cos $2 (Ig ® Ox/2) + sin $2 (Ig ® oy/2)| 
+2rJo,@o,/4, (2) 


where w1; and ġ; are the control parameters denoting 
the amplitude and the angle in the ry-plane of the exter- 
nal rf fields, respectively. Above J denotes the spin-spin 
coupling constant. Typically we have wi; > J, which 
justifies the above assumption of negligible one-qubit op- 
eration time compared to two-qubit operation time. Ne- 
glecting one-qubit operation time in evaluating the exe- 
cution time amounts to identifying the matrices U; and 
Uz which differ by an element of K = SU(2) & SU(2). 
Thus the relevant space for evaluating the time opti- 
mal path is the coset space SU(4)/SU(2) ®@ SU(2). To 
find the time optimal path connecting the unit matrix I 
and the matrix Uaig, therefore, amounts to finding the 
time optimal path connecting cosets [I] and [Uaig], where 
[U] = {kU|k € K}. The Lie algebra su(4) of SU(4) is 
decomposed as su(4) =€@p uA, where 


£ Span({il ® o;/2,i0; ®I/2}), (j = x,y,z) (3) 
p = H= Span({ia; Q on/4}), (j,k =2,y,2). (4) 


Note that they satisfy the commutation relations 


[p,p] C È. (5) 


The decomposition of a Lie algebra g into € and p, sat- 
isfying the above commutation relation, is called the 
Cartan decomposition. The Cartan subalgebra b = 
Span({ia; ® o;/4}) C p plays an important role in the 
following construction. A general theorem of Lie algebras 
proves that any element Uag E€ SU(4) has a KP decom- 
position Uag = kp,k € K = expt and p € P = expp. 
Moreover, any matrix p € P is rewritten in the conjugate 
form p = kihki, where kı € K and A is an element of 
the Cartan subgroup H of SU(4), 


ee C £, [p, €] C p, 


H = exph = 4 exp [7 Da 0; ® 0; 
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(6) 
Therefore we have a corresponding Cartan decomposi- 
tion for a group element as Uag = k2hkı, where k; € K 


and h € H. The implementation of a quantum algo- 
rithm based on a Cartan decomposition requires shorter 
execution time in general d id Gal. 

To be more concrete, let us consider implementing two- 
qubit Grover’s database search algorithm U;; with an 
NMR quantum computer. The data is encoded in one 
of the basis vectors |00),|01),|10),|11) and the gate U;, 
picks out a particular binary basis vector |ij) as a “target 
file” . Here we restrict ourselves within U10 which 
picks out the file |10) with a single step. The unitary 
matrix representing this algorithm takes the form 


cto p 
0 0 i 

MO) 4 gp? gs (7) 
0 O-1 0 


The Cartan decomposition of an arbitrary U € SU(4) is 
carried out_explicitly as follows. We first introduce the 
Bell basis 


|Yo) = (1/v2)(|00) + |11)), 
|W) = (i/V2)(|01) + |10)), (8) 
|W2) = (1/v2)(|01) — |10)), 
|W) = (i/v2)(100) — |11)). 


The transformation rule of a matrix U with respect to the 
standard binary basis |00),|01),|10),|11) into that with 
the Bell basis |W;) is U + Ug = QUQ, where 


(9) 


The matrix Q defines an isomorphism between K = 
SU(2) @ SU(2) and SO(4) and is used to classify two- 
qubit gates Ep Namely, it is easy to verify that 
QikQ € SO(4) for k € K. Moreover, Q diagonal- 
izes the elements of the Cartan subgroup, viz QRQ = 
diag(e%°, #1, e12, e183) for h € H. Therefore we find for 
U= ko hkı that 


Up = QUQ = QtkoQ - QthQ - Qtk2Q = OzhpO1, 


where O; = Q'k;Q € SO(4) and hp = QthQ is a diag- 
onal matrix. From UZUg = OTh? O1, we notice that 
UŁUp is diagonalized by O; and its eigenvalues form 
the diagonal elements of h. Finally O2 is found as 
Oz = Up(hpO,) "1. 

We apply the above strategy to find the Cartan decom- 
position Ujg = kəhkı. An example of the time optimal 
control is 


ky = 1@l, 
h = el(t/4) (02802 -0yB0y) (10) 
ka = e7 r/4)oz Q eilt/2V2)(o2+oy) 


To implement this decomposition with an NMR quantum 
computer, such terms as e(*/4)(%#®¢=) must be rewritten 


TABLE I: time optimal pulse sequences for Grover’s algorithm Uio and the warp-driven algorithm, W4Ui9o. The hydrogen 
nucleus is the qubit 1 while the carbon nucleus is the qubit 2. Here X (Xm) and Y (Ym) denote 7/2-pulse around x (—2) and 
y (—y) axis, respectively. The symbol Pi(@) denotes a 7-pulse around a vector (cos 0, sin 9,0) in the Bloch sphere. The symbol 
(1/2) indicates the length of the idle time, during which no external pulses are applied. The number of pulses is reduced from 


10 to 4 and the execution time is halved by adding the extra gate W4. 


Gate Pulse sequence 
Uio 1: X (1/27) Xm Y (1/2J 
2: X (1/2J) Xm Ym 1/2J 
WaU io 1: Xm Pi(—7/4) 
2: 


in favor of the subset of generators of SU(4) in the Hamil- 
tonian (1). We verify, for example, that 


eiT (62802) = [eile g eiT Doy eilt 02802) 


x [er /A) oe Q ell 80u], (11) 


The extra gate W which possibly shortens the execu- 
tion time must be simple enough so that we can deduce 
the output of Uag efficiently from that of WU,), using 
classical computation only. We call such gates warp-drive 
gates. Since a matrix Uag is an element of a compact 
group U(2”), there always exist such warp-drive gates 
which will reduce the execution time. In the present 
work, we choose the permutation matrices of the binary 
basis vectors as candidates for such gates. Since there 
are four basis vectors, there are 4! = 24 permutation ma- 
trices. Note that we are not required to examine all of 24 
permutations since 18 of them are obtained by applying 
one-qubit rotations on the following six permutations: 


12 21 
Wo = I4, Wi = Ucnor, W2 = Uenor: 
— y2 21 
W3 = Uswap; W, = Ucp = UcnotUénot: (12) 
_ 772 — 7721 12 
Ws = Us = UcnorUecNnor: 


where Uswap represents the swap gate and Uġnor is 
the CNOT gate whose control bit is ¿ while the target 
bit is j. The matrices Usp and UZ, have been utilized 
to cyclically permute the state populations to generate 
pseudopure states 19}. 

We have optimized the execution time of the matrices 
WiUi0 (0 < i < 5) by utilizing the Cartan decomposition 
outlined above and found that the execution time is 1/J 
for i = 0,1 and 2 while 1/2J for i = 3,4 and 5. A time 
optimal control for the warp-driven gate W4Uj9 = kehk, 
is 


ky = etl /2V2)(o2—22) Q ei(t/4)oe 

h = ere ree, (13) 

ko = IQ grea! 
Therefore the execution time T satisfies 7JT/2 = 7/4, 
yielding T = 1/2. Since the decomposition (13) con- 
tains generators which do not exist in the Hamiltonian, 


a trick similar to (I) must be employed. The one-qubit 
gates are realized using rf-pulses in an NMR quantum 


) 
) 
1/23) 
1/2J) Pi(m) 


Execution Time 
X Ym 1/J 
Y Pi(m/4) 
X 1/2J 


computer. Table I shows the actual NMR pulse se- 
quences derived from the Cartan decompositions of U10 
and W4Uj9. We call the hydrogen nucleus and the car- 
bon nucleus as qubit 1 and qubit 2, respectively, in Table 
I and the rest of this Letter. In spite of the extra gate 
W4, the latter requires less pulses and half of the execu- 
tion time required for the former. The permutation W4 
maps the binary basis vectors as W4 : |10) ++ |11), and 
accordingly the output state of W4U10 upon acting the 
initial state |00} is |11). 

In our experiments, a 0.6 ml, 200 mM sample of 
carbon-13 labeled chloroform (Cambridge Isotope) in d-6 
acetone has been employed as a two-qubit molecule and 
data were taken at room temperature with a JEOL ECA- 
500 NMR spectrometer, whose hydrogen Larmor fre- 
quency is approximately 500 MHz. The measured spin- 
spin coupling constant is J = 215.5 Hz and the transverse 
relaxation time is Tə ~ 7.5 s for the hydrogen nucleus and 
Tə ~ 0.30 s for the carbon nucleus. The longitudinal re- 
laxation time is measured to be Tı ~ 20 s for both nuclei. 
The initial state |00) is prepared as a pseudopure state 
generated by the field gradient method Ro. 

Figure 2 shows our experimental result. The spectra 
for the 13C nucleus (qubit 2) are shown in the panels. The 
spin states of both nuclei are inferred from these spectra. 
The position of the peak depends on the state of the qubit 


100) 


` 77.49 ppm 
111) 


` 79.20 ppm 


FIG. 2: NMR spectra obtained by applying a reading pulse to 
the °C nucleus. The rear row shows the pseudopure initial 
state |00) while the middle row indicates the state |10} ob- 
tained by executing U10 on |00). Finally the front row shows 
the spectrum obtained from warp-drive quantum computing 
Wa4U10, where the peak corresponds to |11). 


1: peaks at 77.49 ppm (part per million) and 79.20 ppm 
stand for the states |1) and |0), respectively, of the qubit 
1. The positive (negative) amplitude indicates qubit 2 
being in the state |0) (|1)). Thus the peak in the rear row 
shows the initial pseudopure state is |00) while the peak 
in the middle row tells us that the execution of Grover’s 
algorithm U10 generates the state |10). Finally the front 
row shows the spectrum obtained after the execution of 
the warp-driven gate W4U10. The peak shows that the 
resulting state is |11), from which the output of U10100} is 
easily deduced as |10). The spectrum obtained with the 
warp-driven gate is sharper than that with the original 
Uio gate. The undesirable peak at |00) is also improved 
by the warp-drive. Smaller number of gates and reduced 
execution time account for these improvements. 

We have also examined the effect of the warp-drives 
for Grover’s algorithms other than Uj9 and found that 
the execution times remain unchanged for Wo,1,2 while 
they are halved for W345. It should not be expected, 
however, that the latter gates reduce the execution time 
of an arbitrary two-qubit gate: clearly, such a “universal 
warp-drive” which reduces execution time of an arbitrary 
quantum algorithm does not exist. In view of the com- 
pactness of the unitary group, however, there exists a fi- 
nite set of simple unitary transformations, such that the 
execution time of any quantum algorithm is reduced by 
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In summary, we have demonstrated both theoretically 
and experimentally that a quantum algorithm may be ac- 
celerated by adding an extra gate to it. We took advan- 
tage of the compactness of the group SU(4) and found the 
time optimal control parameters utilizing the Cartan de- 
composition for it. When applied to two-qubit Grover’s 
algorithm in the NMR experiment the execution time is 
found to be halved along with a reduction in the number 
of the gate pulses from 10 to 4. As a result, application of 
warp-drive quantum computing sharpens the NMR spec- 
trum and reduces the spurious peak, implying reduction 
in decoherence and gate operation errors. How reduc- 
tions in the execution time and the number of gate pulses 
scales with the number of qubits remains a challenging 
problem. 
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